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ABSTRACT 

We present our survey for optically faint variable objects using multi-epoch (8 — 10 epochs over 2 — 4 
years) j'-band imaging data obtained with Subaru Suprime-Cam over 0.918 deg^ in the Subaru/XMM- 
Newton Deep Field (SXDF). We found 1040 optically variable objects by image subtraction for all 
the combinations of images at different epochs. This is the first statistical sample of variable objects 
at depths achieved with 8-lOm class telescopes or HST. The detection limit for variable components 
is ~ 25.5 mag. These variable objects were classified into variable stars, supernovae (SNe), and 
active galactic nuclei (AGN), based on the optical morphologies, magnitudes, colors, and optical- 
mid-infrared colors of the host objects, spatial offsets of variable components from the host objects, 
and light curves. Detection completeness was examined by simulating light curves for periodic and 
irregular variability. We detected optical variability for 36 ± 2% (51 ± 3% for a bright sample with 
i' < 24.4 mag) of X-ray sources in the field. Number densities of variable obejcts as functions of time 
intervals At and variable component magnitudes i'^^^.^ are obtained. Number densities of variable stars, 
SNe, and AGN are 120, 489, and 579 objects deg~^, respectively. Bimodal distributions of variable 
stars in the color- magnitude diagrams indicate that the variable star sample consists of bright {V ~ 22 
mag) blue variable stars of the halo population and faint {V ^ 23.5 mag) red variable stars of the disk 
population. There are a few candidates of RR Lyrae providing a possible number density of ^ 10^^ 
kpc~'^ at a distance of > 150 kpc from the Galactic center. 

Subject headings: stars: variables: other — supernova: general — galaxies: active — surveys 
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1. INTRODUCTION 

It is well-known that there are a wide variety of objects 
showing optical variability by various mechanisms in var- 
ious time scales in the universe. For example, stars in the 
instability strip in a color-magnitude diagram show peri- 
odic variability by pulsations (Cepheids and RR Lyrae), 
while cataclysmic variables show emergent variability in- 
duced by accretion flows from companion stars (novae), 
magnetic reconnection or rotations (dwarf stars). Ther- 
monuclear or core-collapse runaway explosions of massive 
stars are observed as supernovae (SNe). Recent stud- 
ies on optical afterglows of gamma-ray bursts (GRBs), 
which are extremely energetic phenomena, are starting 
to reveal their nature. Active galactic nuclei (AGN) 
show irregular flux variability not only in the optical 
but also in other wavelengths. Optically variable ob- 
jects described above change their intrinsic brightness, 
but observed variability can be caused by other reasons. 
Searches for microlensing events have aimed to test the 
hypothesis that a signiflcant fraction of the dark matter 
in the halo of our Galaxy could be made up of MAs- 
sive Compact Halo Objects (MACHOs) such as brown 
dwarfs or planets. Many asteroids in solar system have 
been found as moving objects. Many surveys for these 
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optically variable objects have been carried out and have 
contributed to many important astronomical topics. 

The cosmic distance ladder^ for example, were estab- 
lished partly based on studies of optically variable ob- 
jects which are standardizable candles. Observations 
of variable stars in nearby galaxies enable us to deter- 
mine distanc es to them by period-lumin osity relations 
of Cepheids (jSandage fc Tammannll2006l . and references 
therein). Recent deep high-resolution monitoring obser- 
vations with Hubble Space Telescope (HST) increased 
the number of nearby galaxies whose distances were mea- 
sured using Cepheids. Well-sampled light curves of type 
la supernovae (SNe la) in galaxies with distance mea- 
surements by Cepheids (jSaha et al.| [2006) showed empir- 
ical relationships between their light curve shapes and 
their luminosity: intrinsica lly brighter SN e la declined 
more slowly in brightness (jPhillipgl I1993D . These rela- 
tions led to the discoveries of an accelerated expansion 
of the universe by two i ndependent SN la search teams 
(jPerlmutter et al.l Il998t iRiess et all 119981 ) . In ongoing 
and planned SN surveys, SNe la are expected to set tight 
constraints on the value of the dark energy and its cos- 
mological evolution if any. 

Variable stars have been also used to trace Galac- 
tic structure as standard candles. RR Lyrae are old 
(> 8 Gyr) bright {My ~ 0.6 mag) low-mass pulsat- 
ing stars with blue colors oi B — V ^ 0.3 and are ap- 
propriate for tracing structures in the halo. They can 
be selected as stars showing large variability (0.5 — 1.0 
mag in F-band) in short periods of 0.3 — 0.5 days 
(jVivas et al.ll2004D . RR Lyrae selected using multi-epoch 
imaging dat a obtained by th e Sloan Digital Sky Sur- 
vey (SDSS; lYork et all l2000l ) provided a possible cut- 
off of the Galactic halo at_~ 65 kpc from the Galac- 
tic center ([Ivezic et al.l |2000[ ) . The sub sequent stud- 
ies on RR Lyrae by SDSS (jlvezic et al . 2004a, 200J; 
ISesar et aD l2007l) and the QUasa r Equa torical Survey 
Team (QUEST: iVivas et al.l l200l rvTvai^ fc Zinm l2006fl 
found many substructures of the Galactic halo includ- 
ing those already known. 

Optical variability of t he first luminous AGN (quasar) 
3C 273 were recog nized (ISmith fc HofHeitl ll963D just af- 
ter its discovery ([Schmidt! Il963f l. Since then, it has 
been known that almost all quasars show optical vari- 
ability; indeed, it is one of common characteristics of 
AGN. Quasar surveys such as SDSS and the 2-degree 
Field Quasar Redshift Survey (20Z: iBovle et al.ll2000[ ) 
have been carried out mainly using optical multi-color 
selection s. These surveys have made large catalo gs of 
quasars ()Croom et al.l l200l iSchneider et all l2005f ) and 
also found very high-z quas ars close to the reioniza- 
tion epoch (F an et aI1l2006al lbl). But efficiency of find- 
ing lower-luminosity AGN using color selections are ex- 
pected to be low because AGN components get fainter 
compared with their host galaxy components. On the 
other hand, recent deep X-ray observations with ASCA, 
XMM-Newton, and Chandra satellites effectively found 
many distant low-luminosity AGN as well as obscured 
AGN. This high efficiency owes to faintness of host galax- 
ies in X-ray and high transparency of dust to X-ray pho- 
tons. X-ray observations can identify faint AGN easily 
while finding even unobscured AGN is difficult using op- 
tical color selections. However, since deep X-ray surveys 
over wide fields require a lot of telescope time (e.g. 2 Ms 



expo sure in the Chandra Deep Field-North: iBrandt et all 
1200 It ), optical variability is being recognized again as a 
good tracer for AGN. Some studies ha ve succeeded in de- 
tecting optica.1 varia.bility of quasa rs ()Hook et al.lll994l : 
iGiveon et aP 11999': 'Hawkins '2002), and finding many 
quasars by optical variability with th e completeness as 
high as classical UV exc ess selections (|Hawkins fc Veronl 
[1991 llvezic et al.l[200l . Comparison o f SDSS imag- 
ing data w i th older plate im aging data ([de Vries et al.l 
[20031 [20051 ISesar et al.l [200l an d SDSS spectrophoto- 
metric data ( Vanden Berk et al] |2004| ) showed a clear 
anti-correlation between quasar luminosity and optical 
variability amplitude, as indicated in previous studies 
([Hook et al.lll994l : iGiveon et al.lll999D . Therefore, op- 
tical variability can be an efficient tool to find low- 
luminosity AGN if variable components can be extracted. 
Ultra-deep optical variability surveys with Wide-Field 
Planetary Camera 2 (WFPC2) installed on HST ac- 
tually found 24 galaxies with variable nuclei down to 
Vnuc = 27.5 mag and I^uc = 27 mag, which are as 
faint as ne arby Seyfert galaxies ( —15 < Mb < —19) 
at z ~ 1 ([Saraiedini et all l2000l [2001 . ICohen et al.l 
([2006( 1 also found several tens of AGN by optical vari- 
ability using multi-epoch data with Advanced Camera 
for Survey (ACS) on board HST. The number densi- 
ties of AGN selected via optical variability can be of the 
same order as those selecte d via deep X-ray observations 
([Brandt fc Hasingeij|2005( ). However, samples from the 
HST imaging data were not large enough for statistical 
studies. The comparable number densities of variability- 
selected AGN may indicate that there could be several 
populations of AGN with different properties. Furthre- 
more, a search for faint tra nsient objects with Suprime- 
Cam ([Mivazaki et al.l[2"002l ) on Subaru telescope revealed 
very faint AGN variability in the nuclei of apparently 
normal galaxies, and found that such nuclei show very 
rapid (~ a few day s) nuclear variab ility with a large 
fractional amplitude ([Totani et al.ll2005i) . Such behavior 
is similar to that of the Galactic center black hole, Sgr 
A*, rather than bright AGN, possibly indicating different 
physical nature of accretion disks of very low luminosity 
AGN and very luminous AGN. 

Object variability could affect searches for ra re objects 
using non-si multaneou s obser vational data, live et al.] 
([20060 and lOta et all ([2007D searched for Lyman-a 
emitters (LAEs) at z ^ 7.0 by comparing narrow- 
band data with broad-band data obtained more than 
a few years earlier. They succeeded in identifying a 
bright candidate spectroscopically as a LAE at z = 
6.94. Another candidate without spectroscopic identi- 
fications could be just a transient object and the au- 
thors treated it as a marginal can didate. SN rate 
studies based on optical v ariabilit y (iPain et aTl 120021: 
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20061: iPoznanski et al.ii2007i : lOda et al.ii2007i ) have some- 



times confronted problems in classifying events due to in- 
sufficient observational data such as spectroscopic identi- 
fications, sufhcient time samplings, and multi-wavelength 
measurements. Given limited time on large telescopes, 
one cannot make spectroscopic observations of all vari- 
able candidates; thus, multi-wavelength imaging data 
and light curves in year-scale baselines are very useful 
to separate SNe from other kinds of variable objects. 
Recently, many optical variability surveys with dense 
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time samplings have been conducted for various astro- 
nomical purposes. Dividing exposure time into mul- 
tiple epochs in extremely deep surveys using large or 
space telescopes also enables us to explore optical faint 
variability. Projects conducted recently or ongoing 
are Faint Sky Variability Survey fFS VS: iGroot et all 
I2003D . Deep Lens Survey (DLS; B ecker et all ' 



Supernova Legacy Survey (SNLS; lAstier et al 
Great Obserya,tory Origins Deep Survey (GOODS; 
IStrolgcr et all I2004D. Grot h-Westphal Survey Strip 
(GSS; Saraiedi ni etall [200 61) . Hu bble D eep Field-North 
(HDF-N; Saraiedini et al. 200 0. j2003l) Hubble Ultra 
Deep Field (HUDF; ,Cohen et al]l2006f ). Sl oan Digital 
Sky Survey-H (SDSS-H) Supernova Survey (|Sako et all 
l2007f) . Densely sampled observations targeting GRB 
orphan afterglows also have been carried out (e.g., 
iRau. Greiner. fc Schwar j I2006D . We focus here on 
our survey for optically faint variable objects since 
2002 in the Subaru/XMM- Newton Deep Field (SXDF, 
iSekiguchi et all I2004I . [20071 ). The data have been taken 
by the Subaru/XMM-Newton Deep Survey (SXDS) 
project. Our study described in this paper provides the 
first statistical sample of optically faint variable objects 
and is unique among studies using 8-lOm class telescopes 
and HST in its combination of wide field coverage and 
depth. The Suprime-Cam imaging data and the anal- 
ysis for finding optically variable objects are described 
in [21 followed by other observational data in SJH] We 
describe object classifications in fJH and detection com- 
pleteness in 5j5l The results are discussed for statistics 
of the whole sample in f^Hl and for variable stars alone in 
[JT] We summarize our results in Sj8l In this paper, we 
adopt throughout the AB magnitude system for optical 
and mid-infrared photometry. 

2. OPTICALLY FAINT VARIABLE OBJECT SURVEY WITH 
SUBARU SUPRIME-CAM 

2.L Subaru/XMM-Newton Deep Field (SXDF) 

We are carrying out an optically faint variable object 
survey using a wide-field camera, Suprime-Cam, on the 
prime focus of Subaru 8.2-m telescope. The widest field- 
of-view (34' x 27') among the optical imaging instru- 
ments installed on HST and 8-IOm telescopes provides 
us a unique opportunity for statistical studies of rare and 
faint objects which are less affected by cosmic variance 
than those derived from smaller fields. 

The SXDF is centered on (02'^I8'"00% -05:00:00) 
in J2000.0, in the direction towards the Galactic 
halo, (/,6) = (169°, -60°). The SXDF is an on- 
going multi-wavelength project from X-ray to radio 
JSckiguchi et al. 2004, 2007), exploring aspects of the 
distant universe such as the nature of the extragalac- 
tic X-rav populations (l\yatson et all 120051 : lUeda et all 
120071: lAkiyama et al.ll2007D, la rge-scale structures at high 
redshift dOuchi et al.' '2005a"b') and cos mic history of 
mass assembly (Kodania ct al. 2004; Yani ada et al.|[2005l : 
iSimpson et al][2006bi) . The field coverage is ~ 1.2 deg , 
which consists of five pointings of Suprime-Cam (SXDF- 
C, SXDF-N, SXDF-S, SXDF-E, and SXDF-W, respec- 
tively, see Figure [T]). We use the i'-band-selected photo- 
metric catalog s in B, V, Rg, i ' and z'-bands, which are 
summarized in lFurusawa et al.l ()2007). The depths of the 
catalogs which we use are almost the same among the five 
fields, 28.4, 27.8, 27.7, 27.7, and 26.6 mag (3cr, 2'.'0(f>) in 



B, V, Rc, i' and z'-bands, respectively (Table [T]). Op- 
tical photometric information used in this paper except 
for variability measurements is derived from these cat- 
alogs and are not corrected for the Galactic extinction 
in the SXDF d irection, E(B - V) = 0.019 - .023 mag 
(jSchlegel et al.lll998. see lFurusawa et al.ll2007[ ). 

2.2. Subaru Suprime-Cam Imaging Data 

Our survey for optically faint variable objects in 
the SXDF is based on multi-epoch i'-band Suprime- 
Cam i maging data, starting in S eptem ber 2002 
(jFuri^awa et aQ 120071 : lYasuda et all l2007f ). The 
Suprime-Cam observations are summarized in Table |2l 
These observations include those for an extensive high- 
z SN search using well-sa mpled i'-ban d images i n 2002 
(|Yasuda et al.ll2003L l2007[ l with which iDoi et all (|2003l ) 
reported the discoveries of 13 high-z SN e. In addition to 
observ ations used to make the catalogs in lFurusawa et ah! 
(|2007| ). we carried out i'-band observations in 2005. The 
final images were taken in October 2003 in the two fields 
(SXDF-N and SXDF-W) and in September 2005 in the 
remaining three fields (SXDF-C, SXDF-N, and SXDF- 
E) . The numbers of the observational epochs (iVopoch) in 
the five fields are 10, 8, 8, 10, and 8, respectively. Time 
intervals of the observations are from 1 day to 3 years in 
observed frame. 

In each epoch, stacked images were made in a stan- 
dard method fo r the Supr ime-Cam data using the NEKO 
software (.Yagi et all l2002i) and the SDFRED package 
(lOuchi et al.l I2004D . Some images taken on different 
dates were combined together to make all the depths 
almost the same. Then, we geometrically transformed 
the images using geomap and geotran tasks in IRAF^^ in 
order to match the coordinates in better accuracy each 
other. The root-mean-square of residuals of the coordi- 
nate differences are typically 0.4 pixel. 

Typical exposure times of one hour provide limiting 
magnitudes of = 25.2 — 26.8 mag (< >mcdian= 
26.0 mag, 5(7, 2'.'0(/)) . The fuU-width-at-half-maxima 
(FWHM) of point spread function (PSF) were 9 = 0'.'52- 
1'.'54 (< 6 >mcdia„= (^.'70, see Tabled]) 

For our variability study, we used regions overlapping 
in all the epochs in each field. Considering saturation and 
non-linearity of CCDs at high signal levels, we excluded 
regions around bright objects to ensure reliable detection 
of object variability. The total effective area is 0.918 deg^ 
(Table [1]) shown as light gray and dark gray regions in 
Figure [TJ 

2.3. Variability Detection 

The Suprime-Cam images have different PSF size (and 
shape) in every epoch because of time-varying atmo- 
spheric seeing. Therefore, we can not measure variability 
of objects b y simply compar i ng flu x wit hin a f ixed small 
apertu re as lSaraiedini et all (|2000[ ) and ISaraied ini et aO 
$mf) did for the HST WFPC2 images. ICohen et all 
( 2006) used total magnitude for variability detections in 
the images obtained with the Advanced Camera for Sur- 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

Pixel scale of Suprime-Cam is 0''202 pixel^^. 
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relative RA [orcmln] 

Fig. 1. — The Suprime-Cam imaging fields (light gray and dark 
gray regions) in the SXDF, superimposed XMM-Newton EPIC 
imaging fields (black solid line) and Spitzer IRAC imaging fields 
(light gray region). The coordinates are measured relative to the 
center of the SXDF, (02hl8™00^ -05:00:00) in J2000.0. Our op- 
tically variable object survey has been carried out in the entire 
Suprime-Cam fields over 0.918 deg^ except for regions around 
bright objects (white blanks). The five pointings of the Suprime- 
Cam, SXDF-C, SXDF-N, SXDF-S, SXDF-E, and SXDF-W, are 
described. 
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Fig. 2. — Examples of variability detection. In each row, an im- 
age on a reference date, an image on another date, and subtracted 
image are shown from left to right. Variable objects are clearly 
seen as positive or negative residuals in the right panels. All these 
three objects were spectroscopically identified. The top object is a 
M dwarf star, the middle object is a broad-line AGN at z = 2.150, 
and the bottom object is a SN la at 2 = 0.606. The spectra and 
light curves are shown in the top row of Figure \3\ the third row of 
Figure |4l and the bottom row of Figure [5] respectively. The box 
sizes are 20". 



vey (ACS) on HST because the PSF varied with loca- 
tions on the CCDs of ACS. Accurate photometry itself 
of any kind for extended objects is often difficult. Since 
variable components are almost always point sources, 
the light from non-variable extended components just 
increases the noise for variability detections. Hence, 
we used an ima ge sub traction method introduc ed by 
lAlard fc LuptonI (|1998D and developed by lAlardI (|200Clf ) 
instead of using total magnitudes of objects. This sub- 
traction method using space- varying convolution kernels, 
which are obtained by fitting kernel solutions in small 
sub-areas, enables us to match one image against an- 
other image with a diff'erent PSF. We can then detect 
and measure variable objects in the subtracted images. 
We applied this method for all the possible two pairs 
of the stacked images at different epochs and detected 
variability of the objects in each field. Figure [2] shows ex- 
amples of the image subtration for three of the variable 
objects. 

We now describe details of the procedures for vari- 
ability detection and our definitions of variable objects. 
First, we smoothed all the subtracted images with the 
PSFs and detected positive and negative peaks of surface 
brightness in subpixel unit. Then, in all the subtracted 
images, we selected objects whose aperture flux in a fixed 
diameter of 2'.'0 were above Stibf or below — 5crbf of the 
background fluctuations (cTbf) within the same aperture. 
The background fluctuations Ubf is obtained by calculat- 
ing the standard deviation of 2'.'0 aperture flux in ~ 200 
places around each positive or negative peak in the sub- 
tracted images. We inspected these variable candidates 
visually to exclude candidates around mis-subtracted re- 
gions. Second, we deflned the central coordinates of the 

The number of the combinations is m , C2 . 



variable object candidates as the flux-weighted averages 
of all the coordinates where significant variability were 
detected in the subtracted images in the first procedure. 
Third, we did 2'.'0-diameter aperture photometry for all 
these variable object candidates in the subtracted im- 
ages, the reference images of which were an image with 
the smallest PSF size in each field^°, and scaled the aper- 
ture measurements to the total flux assuming that the 
variable components were always point sources. The fac- 
tors to scale 2'.'0 aperture flux to total flux were calcu- 
lated as ratios between 2'.'0 aperture flux and total flux 
for point sources in each image (typically 1.1 — 1.3). 
We constructed total flux light curves of all the candi- 
dates and deflned optically variable objects as objects 
showing more than 5aic {<7ic'. photometric errors in the 
light curves) variability in at least one combination of 
the epochs in the light curves. We note that some vari- 
able objects have unreliable photometric points in the 
light curves. In the subtraction method, very slight off- 
sets (even ~ 0.05 pixel) for bright objects between two 
images which are used for the subtraction could provide 
signiflcant dipole subtractions. When an object has both 
a positive peak above ba^f and a negative peak below 
—5ahf locally (around ~ 10 pixel square) in a certain 
subtracted image used to make the light curve, the pho- 
tometric point of the object at that epoch are regarded 
as unreliable. These epochs were not used in evaluating 
variability. 

We followed these procedures for all the flve flelds. Fi- 
nally, we made an optically variable object catalog con- 
sisting of 1184 objects in the area of 0.918 deg^ (1290 
variable objects deg~^). However, for widely spread 

20 They are images obtained on 02/09/29 and 02/09/30 for all 
the fields. 
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galaxies or bright objects, subtracted images are some- 
times insecurely extended and show bumpy features al- 
though saturated or non-linearity pixels were removed 
from the detections in advance. Variability of these ob- 
jects are considered marginal. When we exclude such 
marginally variable objects, the number of variable ob- 
jects is 1040 (1133 variable objects deg-^). These 1040 
variable objects have 9236 photometric points in their 
light curves, and 359 unrcliablly measured photometric 
points (3.9%) of 169 objects were not used in evaluat- 
ing variability. In the following discussions, we use this 
secure sample of 1040 variable objects. We note that 
detection completeness is discussed in iJSJ 

In principle, variability should be measured as changes 
of luminosity of objects. To do so, it is necessary to 
measure total flux of variable objects. However, this is 
possible only for point sources such as variable stars and 
quasars whose luminosity is much larger than that of 
their host galaxies. To measure variable AGN separately 
from their host galaxies is difficult using ground-based 
imaging data especially when the flux of variable object 
is faint compared to the flux of the host. We measure the 
flux of the variable component (A/ = /i — /2) as describe 
above. In this paper, we use variable component magni- 
tudes (i^jjj.; = — 2.5 log(|A/|) -I- mzoro) to describe object 
variability, and we adopt the total magnitudes which are 
measured in the SXDF catalogs as the magnitudes of 
host objects. 

2.4. Assignment of Host Objects 

In order to investigate how far the variable components 
are located from the centers of possible host objects, and 
their properties such as magnitudes and colors in the cat- 
alogs, we assigned host objects to all the variable objects. 
In general, AGN and SNe have their host galaxies while 
variable stars have no host objects. In this paper, we 
defined variable stars themselves as host objects to treat 
all the variable objects equally. We calculated projected 
distances of the variable components from the surround- 
ing objects in the stacked images taken when the variable 
objects were in the faintest phases, and defined the near- 
est objects as the host objects. The positions of the host 
objects were measured as peaks of surface brightness of 
the images smoothed by the PSFs. In the faintest phases 
of variable objects which may disappear (i.e. SNe), host 
objects should be detected, but high- 2; SNe which ocuur 
in diffuse galaxies seem hostless in some cases. Assigned 
host objects for such SNe may not be real host galaxies. 
It is also possible that transient objects on faint host ob- 
jects may have no equivalent host objects in the SXDF 
catalogs because the SXDF i'-band images for the cata- 
logs were made by stacking all the exposures except for 
those in 2005. Among 1040 variable objects, the number 
of such objects, which do not have equivalent objects in 
the SXDF catalogs, is 8. These possible transient objects 
may be hostless SNe and are treated as SNe if their light 
curves satisfy a criterion described in i )4.2.2l The number 
of the possible transient objects is small and we do not 
exclude or discriminate them in most of the following dis- 
cussions. In [JSl we discuss number densities of transient 
objects, which are defined as objects not detected in their 
faint phases; objects appear in some epochs at positions 
where there are no objects detected in other epochs. In 
this paper, the number of detections for each transient 



object is not considered in the definition of transient ob- 
jects. 

Information on host objects, such as total magnitudes 
and colors, is derived from integrated images which were 
obtai ned in var i ous e pochs in each band as shown in 
iFurusawa et al.l ()2007[ ) and they are averaged over the 
observational time spans. Therefore, properties of ob- 
jects are appropriate for statistical studies and could be 
inaccurate if we focus on a certain object. 

3. ADDITIONAL OBSERVATIONAL DATA 

The SXDF project has carried out multi-wavelength 
observations from X-ray to radio and follow-up optical 
spectroscopy. In this paper, optical spectroscopic results 
and X-ray imaging data are used to confirm the valid- 
ity of our object classifications. X-ray data is also used 
for evaluating the completeness of AGN variability de- 
tection. We use mid-infrared imaging data to classify 
optically variable objects in S|4l In this section, we sum- 
marize these observations and cross-identifications of op- 
tically variable objects. The summary of the X-ray and 
mid-infrared imaging data is given in Table [TJ 

3.1. Optical Spectroscopy 

Since 2002, the follow-up optical spectroscopic obser- 
vations have been conducted with many telescopes and 
instruments (lYasuda et all 120031: lYa mada ct al. 200l 
Lidman et all 120051: IWatson et al.n 2005i: ,Ouchi et all 
2005allbl : lYasuda et al. 20071 ). Instruments and telescopes 
used are the Two-Degree Field (2dF) on the Anglo- 
Australian Telescope (AAT), the Faint Object Cam- 
era And Spectrograph (FOCAS; Kashik awa et al.ll2002D 
on Subaru tele scope, the Visible M ultiObject Spectro- 
graph fVIMOS: lLe Fevre et al.ll2003l ). and the FOcal Re- 
ducer/low dispersion Spectrograph 2 (F0RS2) on Very 
Large Telescope ( VLT), the Gemini Multi-Object Spec- 
trograph (GMOS; iHook et all [200I on Gemini-North 
telesc ope, and the Echellette Spectrograph and Imager 
(ESI: ISheinis et al.r 2002) on Keck-II telescope. All the 
observations except for the Keck-II ESI, which is a high- 
dispersion echellete spectrograph, have been carried out 
in relatively low spectral resolution mode {R ^ 500). 
Most of the spectroscopic observations were done in 
multi-object spectroscopy mode. 

Out of 1040 variable objects, 119 objects were tar- 
geted for spectroscopy. Of these 119 objects, 99 objects 
were identified and the redshifts were determined. The 
spectroscopic sample includes 3 stars, 14 host galaxies of 
SN e, 6 of which are sp e ctroscopically ident ified as SNe 
la (jLidman et all 120051: lYasuda et al.ll2007D . 53 broad- 
line AGN, 4 [Ne V] emitting galaxies and 25 galaxies. 
Strong [Ne V] emission lines are only observed among 
AGN, which emit strong ionizing photons above 97.2 eV, 
and [Ne V] emitting galaxies are considered to host AGN. 
Other galaxies which appear to be normal galaxies also 
can host not only SNe but also AGN because many AGN 
in our sample are faint compared to the host galaxies 
so that the spectra may not show significant features of 
AGN origin. The redshifts of 96 extragalactic variable 
objects are 0.240 < z < 4.467. We show 11 examples of 
the spectra obtained with FOCAS on Subaru and F0RS2 
on VLT in Figure [3l Figure [H and Figure El The top 
two rows of Figure [3] are those for M dwarf stars showing 
bursts in 2002 with red optical colors oi B — V = 1.34 



6 



T. Morokuma et al. 



and 1.36. The bottom row of Figure [3] is an early-type 
star with B — V = 0.44. We also show spectra of four 
broad-line AGN at z = 0.867, 1.066, 2.150, and 3.553 in 
Figure m and those of SN host galaxies at z = 1.239, 
0.505, 0.517, 0.606 in Figure [H The light curves are also 
shown in the figures. The photometric points in these 
light curves are measured as differential flux in the sub- 
tracted images. 
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Fig. 3. — Spectra in arbitrary unit and light curves of three vari- 
able stars. All the spectra were obtained with FOCAS on Subaru 
telescope. The observational configurations were the 300B grism 
and SY47 order-sort filter for the top two stars and the 150 grism 
and SY47 order-sort filter for the bottom star. The top two rows 
are M dwarf stars with B — V = 1.34 and 1.36. The bottom row 
is an early- type star with B — V = 0.44. The light curves of dif- 
ferential flux were plotted in filled circles in linear scale (in ADU) 
because we can not plot magnitudes in their faintest phases. The 
zeropoint is 34.02 mag. Unreliable photometric points are plotted 
in open circles. 
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Fig. 4. — Spectra in arbitrary unit and light curves of four 
extragalactic variable objects which are classified as AGN based 
on their light curves and variable locations {central variability). 
All the spectra were obtained with FOCAS on Subaru telescope. 
The observational configurations were the 300B grism without any 
order-sort filters for the top and third AGN and the 150 grism 
and SY47 order-sort filter for the second and bottom AGN. The 
redshifts are 0.867, 1.066, 2.150, and 3.553, from top to bottom. 
Dashed lines indicate detected emission and absorption lines. The 
zeropoint the light curves is 34.02 mag. 



Since spectroscopic observations of variable objects 
were mainly targeted for X-ray sources and high-z SN 
candidates, the spectroscopic sample of variable objects 
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Fig. 5. — Spectra in arbitrary unit and light curves of four ex- 
tragalactic variable objects. The top three objects are classified 
as SNe based on their light curves and variable locations {offset 
variability). The bottom object was spectroscopically identified as 
a SN la, SN 2002km (Lidman et al. 200SD, but we can not deter- 
mine whether they are SNe or A GN be cause baselines of their light 
curves arc not long enough (see i|4.2.2l l. The top three objects were 
observed with FOCAS on Subaru telescope. The observational con- 
figurations were the 150 grism without any order-sort filters for the 
top and third objects, and the 300B grism and SY47 order-sort fil- 
ter for the second object. The bottom object were observed with 
FORS2 on VLT using the 3001 grism and t he OG590 order-sort fil- 
ter (see Figure A. 32 in lLidman et al]|2005l for the SN component). 
The redshifts are 1.239, 0.505, 0.517, and 0.606, from top to bot- 
tom. Dashed lines indicate detected emission and absorption lines. 
The zeropoint is 34.02 mag. 



are biased for those classes of variable objects. There- 
fore, we can not discuss much about fractions of variable 
objects only from the spectroscopic sample. 

3.2. X-ray Imaging 

X-ray imaging observations in the SXDF were carried 
out with the European Photon Imaging Camera (EPIC) 
on board XMM-Newton satellite. They consist of one 
deep (~ 100 ks) pointing on the center of the SXDF 
and six shallower (~ 50 ks) pointings at the surround- 
ing regions. In total. X-ray imaging covers most of the 
Suprime-Cam imaging fields (see Figure [1]). The details 
of the XMM-Newton EPIC observations, da ta analyses, 
and o ptical identifications a.re de scribed in lUeda et al.l 
(|2007t ) and lAkivama et all (|2007[ ). The limiting fluxes 
are down to 1 x 10~^^ erg~^ cm~^ s~^ in the soft band 
(0.5-2.0 keV) and 3 x lO^^^ erg"i cm^^ g-i ^^^^ i^^j-^j 
band (2.0-10.0 keV), respectively. 

In this paper, we use a sample which was detected 
with likelihood larger than nine in either the soft or hard 
band and also within our variability survey fields. After 
excluding objects in the regions not used for variability 
detection (Table[T|), we have 481 X-ray sources over 0.808 
deg^ Among 936 variable objects within the X-ray field, 
172 objects were detected in X-rays (165 and 91 objects 
detected in the soft and hard bands, respectively, see 
Table [1]). 

3.3. Mid-Infrared Imaging 

Mid-infrared imaging data in 3.6/im, 4.5/xm, 5.8/xm, 
and 8. 0jini-bands with the InfraRed Array Camera 
(IRAC: lFazio et al.ll2004l ) on Spitzer space telescope were 
obtained in the SXDF as a part of the Spitzer Wide-area 
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Infra Red Extragalactic (SWIRE; iLonsdale et"aI1 12003| . 
[2onl survey. The IRAC data covers almost the entire 
Suprime-Cam field except for a part of SXDF-W. The 
covered field is 0.889 deg^ (97% of the Suprime-Cam 
field) shown as gray regions in Figure [TJ The reduced 
IRAC images were obtained fro m the SWIRE Archiy e 
and the catalogs were made by lAkivama et al.l (|2007D . 
They used the IRAC imaging data for optical identifica- 
tions of X-ray sources. In this paper, we use only the 
3.6/im-band data for the object classification in 52) The 
limiting magnitude is ~ 22.0 mag (total flux, Scr for point 
sources). Among 1040 variable objects, 1028 objects are 
within the IRAC field and 995 objects are detected in the 
3.6/im-band (Table [1]). Spatial resolution of the IRAC 
3.6/im-band is not high, ^ 1", but we can assign IRAC 
identifications to almost all the optically variable objects 
with good accuracy. 

4. OBJECT CLASSIFICATION 

The variable object sample includes several classes of 
variable objects. In this section, we classify these variable 
objects into variable stars, SNe, and AGN using optical 
and mid-infrared imaging parameters. The procedures 
described in i )4.1l and ij4.2l are summarized in a flow chart 
of Figure [H 
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Fi g. 6. — F low chart of object classification. Sec the text in 
i|4.1l and i|4.2l for details. If an object satisfies a criterion at some 
point along the path, it goes to the left; if not, to the right. In the 
first criteria, the IRAC data is necessary, so that we can classify 
variable objects in 97% of the overall field. The last criteria, using 
light curves for extragalactic objects, requires time baselines over 
longer than three years, and therefore only objects in the three 
fields, SXDF-C, SXDF-S, and SXDF-E can be classified. 



We note again that the object photometric information 
is averaged over the observational time spans, but its 
effect on our statistical classifications and discussions is 
small even if objects are variable. 

4.1. Star Selection 

We extracted variable stars from the variable object 
sample using the Suprime-Cam optical imaging and the 
IRAC mid-infrared imaging data. Since the IRAC data 
covers 97% of the variability survey field in the SXDF, 
we concentrate on the 1028 variable objects in the over- 
lapped region of 0.889 deg^ for the following discussions 
such as object classifications and number densities. Star 



selection is based on the optical morphologies, optical 
magnitudes, optical colors, and optical and mid-infrared 
colors of the variable objects. We used total magnitudes 
in both optical and mid-infrared wavelengths to make 
this selection. 

4.1.1. Criteria on Optical Morphology, Magnitude, and 

Color 

First, we assigned to all the variable objects mor- 
phological scores (mgcoro = — 10) and color scores 
(cscore = 0-00 — 1.00) using a method described in 
iRichmondl (|2005f ). The morphological scores niscmc were 
calculated based on magnitude differences S = m2 — rn^ 
between 2'.'0 aperture magnitudes (7712) and 3'.'0 aperture 
magnitudes {m^), and the CLASS_ST AR value s in all the 
five b road-bands from the SExtractor (|Bartin fc Arnout^ 
Il996l ) outputs in the SXDF catalogs. If an object has 
0.10 < 5 < 0.20, we add 1 to m^core- If a CLASS.STAR 
value of an object is larger than some value (typically 
0.90, small variations from field to field) we also add 1 
to mscorc- The color scores Cscorc were assigned based 
on the distances from the stellar loci in two color-color 
planes, B — V versus V — R and V — R versus R — i', us- 
ing the SExtractor output IsophotMag. The stellar loci 
were defined using colors of a small set of bright stars. 
The color score represents the probability that each ob- 
ject is inside the stellar locus in color-color space using a 
Monte Carlo approach. A measure of the degree to which 
the colors of an object (including their uncertainties) fall 
within the stellar locus at the point of closest approach. 
Values range from 0.0 (color is completely outside the 
stellar locus) to 1.0 (color is completely within the lo- 
cus). Objects with larger values for both scores are more 
likely to be stars, and the most likely stars should have 
'Tiscoro ^ 9 and Cgcorc > 0.90. However, faint, red or blue 
stars are likely to have lower values for both the scores 
because of the low signal-to-noise ratio of photometry 
in some bands. In this classification, stars redder than 
G stars were not used for determining the stellar loci. 
The reddest stellar locus is at i? — z' ~ 1.7. Therefore, 
red stars with R — i' > 1.7 can not have high values of 
Cscorc and we do not use Cscoro values for red stars with 
R — i' > 1.7. Galaxies with similar colors can contam- 
inate into the sample. However, they can be excluded 
by an optical-mid-infrared color selection described be- 
low. Another problem is color differences of stars due 
to differences in their metallicity. Bright stars used for 
determining the stellar loci are mainly disk stars, while 
fainter stars are likely to belong to the halo population. 
Color differences in R — i' caused by the population dif- 
ference appear ai R — i' > 1.0, so, faint red stars with 
R — i' > 1.0 tend to have low Cscore values; we decided 
to classify point-like objects with i? — i' > 1.0 as stars. 
Hence, we first assigned stellar likelihood to all the vari- 
able objects through five criteria using only optical pho- 
tometric information, 

a) stars with TOscoic > 9 and Cscore > 0.90 (point-like 
objects with stellar colors), 

b-1) probable stars with 9 > mscoro > 6 and 0.90 > 
Cscoro > 0.60 (less point-like objects with less stellar 
colors), 

b-2) probable stars with mscorc > 9 and R — i' > 1.0 
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(point-like objects with red colors), 

c) possible stars with R — i' > 1.7 (objects with very 
red colors), 

d) non-stellar objects which do not satisfy any criteria 
above. 

Objects satisfying one of the criteria, a), b-1), b-2), or c) 
are considered as stars. 

In order to check the validity of our star selections 
based on only optical imaging parameters, w e used a pop- 
ulatio n synthesis model (Besangon model; iRobin et al.l 
l2003f ). We investigated the expected distributions of 
stars including non-variable stars in the SXDF which 
were selected in the same criteria in a. B — V versus 
V color-magnitude diagram. The distributions contain 
three sequences and two of them, sequences for younger 
disk population and older halo population are reason- 
ably duplicated by the model. Another sequence of faint 
blue point-like objects has almost the same distribution 
as more extended galaxies in this diagram. Then, we 
concluded that these faint blue point-like objects are 
galaxies, and adopted one more criterion, B — V > 
0.08 X y — 1.59 (dashed line in the top panel of Figure 
[7]), to exclude them 

4.1.2. Criteria on Optical- Mid- Infrared Color 

Second, we investigated optical-mid-infrared colors of 
the variable objects in a i? — i' versus R — ms.g^m 
color-color plane. As shown in previous studies on ob- 
ject distributions i n optic al and mid-infrared colo r-color 
planes (Eiscnhard t et all [2004; Rowan-Robinson et al.l 
120051) . stars and galaxies are more distinctly separated in 
optical-mid-infrared color-color planes than in purely op- 
tical color-color planes because of the very different tem- 
peratures of stars and dust. Figure [7] shows color-color 
diagrams of the variable objects for R — i' versus i' — z' in 
the bottom left panel and R~i' versus R — m^ Q^^ in the 
bottom right panel. A galaxy sequence extends widely 
from bottom to top in the blue side oi R — i' colors in 
the R — i' versus R — m^ etim color-color plane while stel- 
lar sequence shows sharp distribution below the line in 
the figure. We defined variable objects satisfying a color 
criterion, R — m^,^^j^^ < 2.6 x {R — i') — 0.5, as stars. 

4.1.3. Results of Star Selection 

Finally, we combined two independent criteria for stars 
and classified the variable objects into four categories; 

1) reliable stars: variable objects which are selected as 
stars a) and have stellar optical-mid-infrared colors 
(78), 

2) probable stars: variable objects which are selected 
as probable stars b-1) or b-2) and have stellar 
optical- mid- infrared colors (24), 

3) possible stars: variable objects which are selected 
as possible stars c) and have stellar optical-mid- 
infrared colors (5), 

4) non-stellar objects: variable objects which do not 
have stellar optical-mid-infrared colors or are cate- 
gorized in non-stellar objects d) (921). 



20 22 24 26 28 
I'— band [mag] 




R-i' [mog] H-t' [mog] 

Fig. 7. — Color-color diagrams of variable objects in versus 
i' — z' (left panel) and R — i' versus R — ms.s/jm (right panel). 
We plot all the variable objects in small gray triangles, probable 
stars in large open circles, and reliable stars in large filled circles. 
A criterion indicated in a dashed line in the top panel, B — V > 
0.08 xV — 1.59, is adopted to exclude point-like blue galaxies. Dot- 
dashed line in the bottom right panel is a line separating stars and 
galaxies, -R — ma.s/jm = 2.6 X (ij — j') — 0.5. A large ellipse indicates 
a region of SPSS quasar (0.1 < z < 5.2) colors, r — i : r — m^sum, 
in IHatziminaoglou et al.l 1120051) and [Richards et al.l (120061 ) . 



All the objects in 1), 2), and 3) satisfy the optical color- 
magnitude critetion {B - V > 0.08 xV- 1.59). The 
numbers of selected variable objects in each selection 
are given between parentheses. In Figure [71 filled circles 
show 1) reliable stars and open circles show 2) probable 
stars and 3) possible stars. 

We have three spectroscopically identified variable 
stars; an early-type star and two M dwarf stars. All these 
stars are securely classified as stars. The early-type star 
and one M dwarf star are classified as 1) reliable stars. 
Another M dwarf star is classified as a 2) probable star 
with a red color oi R — i' = 1.14, a high morphological 
score of mgcorc = 9 (i.e. a point source), and a low color 
score of Cgcorc = 0.34. 

The stellar and galaxy sequences in the R — i' ver- 
sus R — m^ e^n^ color-color diagram join together around 
R — i' ^ 0.2 and R — ma.g/im ~ 0. Objects in this re- 
gion can be contaminated with galaxies. There are sev- 
eral spectroscopically identified galaxies in the region, 
although they are not recognized as stars due to their op- 
tical extended morphologies. Point sources showing opti- 
cal variability can be luminous quasars, but optical-mid- 
infrared colors of quasars are red enough to be discrimi- 
nate d from stars. All the 35 SP SS quasars at 0.2 < z < 
3.7 (jHatziminaoglou et al.l 120051) and 259 SPSS quasars 
at 0.1 < z < 5.2 ([Richards et all I2006f) have redder 
optical-mid-infrared colors (r — TOa.g^im > 0.35, shown 
as a large ellipse) as well as quasars in our sample than 
stars with similar optical colors. We also examined the 
robustness of star selections by comparing the star count 
in T^-band with that predicted by the Besangon model. 
Stars in the SXPF were selected from the catalogs in the 
same criteria as for the variable objects. Star counts in- 
cluding reliable stars, probable stars, and possible stars 
in the SXPF reasonably agree with the model prediction 
in 21 < y < 24. In the bottom right panel of Figure [71 
there can be seen several R — i > 0.6 point sources with 
stellar optical-mid-infrared colors and low color scores 
Cscore On the stellar sequence (gray circles below the dot- 
dashed line). These low color scores Cscoro may derive 
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from non-simultaneous observations and their vari- 
ability. In this way, our classifications based on colors of 
objects could be inaccurate because of non-simultaneous 
observations, but many facts described above indicate 
that variable stars are securely selected through the cri- 
teria with small contaminations and high completeness. 
The reliable stars, probable stars and possible stars, 107 
objects in total (0.889 deg^, 120 objects deg~^), make up 
our sample of variable stars in subsequent discussions. 

4.2. SN/AGN Separation 

After we selected variable stars, we classified non- 
stellar variable objects into SNe and AGN based on two 
parameters from the optical imaging data: the offset of 
the location of variable component from the host object 
and the light curve. 

4.2.1. Criteria on Variable Location 

First, we examined the offsets between the variable 
components and their host objects. Variability should be 
observed at the centers of the host objects for AGN as 
well as stars, while SNe can explode at any position rela- 
tive to the host objects. In order to estimate the errors of 
the locations of variable components for AGN, we exam- 
ined the offsets for optically variable X-ray sources. Most 
of the X-ray sources are considered to be AGN. Some 
of them can be stars emitting X-ray, but for our pur- 
poses, that is not a problem: variable stars also should 
show variability at their centers. Almost all of the X-ray 
detected variable objects 97%) have spatial offsets 
below 1 pixel (0'.'202), and the offsets distribute with a 
scatter of CTogsct, x-ray — 0.29 pixel. We also examined 
offsets for simulated variable objects used in calculations 
of detection completeness in [J5l and found that their off- 
sets from the located positions range with a scatter of 
CoffsGt.sim = 0.51 pixel, slightly larger than that of X-ray 
sources. We set the threshold between central variability 
(variability at the central positions) and ojfset variability 
(variability at the offset positions) to 1.2 pixel, which is 
Scoffsct (using the average of these two error estimates). 

4.2.2. Criteria on Light Curve 

We also used light curves to discriminate SNe from 
AGN. Bright phases where we detect object variability 
should be limited to within one year for SNe for two 
reasons. First, multiple SNe very rarely appear in a 
single galaxy within our observational baselines. Sec- 
ond, any SNe should become fainter than our detection 
limit within one year after the explosions in the observed 
frame. SNe la, the brightest type except for hypernovae 
and some of type-IIn SNe, can be seen at z ~ 1.4 at 
our detection limit (j(,a^j.j 25.5 mag). One year in the 
observed frame corresponds to five months in the rest 
frame at z = 1.4 due to the time dilation by cosmolog- 
ical redshift. Therefore, SNe should fade by more than 
4 ma g below the maximum brightness (e.g., iJha et all 
|2006| ) . This decrease of brightness is a minimum value of 
observed declining magnitudes of SNe and is almost the 
same as the dynamic range of our detection for variabil- 
ity (i^jjj.; = 20.5 — 25.5 mag). Hence, we defined objects 

For example, V-band observations have been mainly car- 
ried out after those in other broad bands were almost finished 
IIFurusawa et al.ll2007l) . 



which are at bright phases for less than one year, and are 
stably faint in the remainder of the observations, as "ob- 
jects with SN-like light curves". In order to apply this 
criterion for the variable objects, it is necessary for light 
curves to cover more than three years at least. In two 
fields (SXDF-N and SXDF-W), our images cover only 
two years, 2002 and 2003; therefore, variable objects in 
these two fields can not be classified as like or unlike SNe. 
Hence, we concentrate on the other three fields (SXDF- 
C, SXDF-S, and SXDF-E) to discuss number densities 
of SNe and AGN based on robust classifications. Among 
693 variable objects in these three fields, 619 objects are 
classified as non-stellar objects. Out of 619 non-stellar 
variable objects, light curves of 27 objects cover only two 
years because of insecure subtractions in some epochs. 
Therefore, we use 592 objects in the following section 
when separating SNe from AGN. The number densities 
derived below can increase by a factor of 1.05(= 619/592) 
to account for the objects lost to poor subtractions. We 
also note that these criteria provide mis-classifications for 
objects with light curves of poor signal-to-noise ratios. 

4.2.3. Results of SN/ACN Separation 

We combined these two criteria; the position offsets 
and the light curves, and applied them for 592 objects 
over 0.566 deg^ to classify variable objects as SNe and 
AGN. Variable objects with significant (> 1.2 pixel) off- 
sets and SN-like light curves (case 1: 158 objects) are 
highly likely to be SNe, while those with small (< 1.2 
pixel) offsets and non-SN-like light curves (case 2: 228 
objects) are highly likely to be AGN. There are no clear 
contaminations in each case in our sample; we have no 
X-ray sources or spectroscopically identified broad-line 
AGN classified as case 1, and also have no SNe with spec- 
troscopic identifications or spectro scopic redshift deter- 
minations of their host galaxies in lYas uda et al.l (j2007f ) 
classified as case 2. 

Objects with small offsets and SN-like light curves 
(case 3: 160 objects) or significant offsets and non-SN- 
likc light curves (case 4: 46 objects) are hard to judge. 
In our time samplings and baselines, AGN can be classi- 
fied as SNe in terms of light curves and objects in case 3 
can be either SNe or AGN. Out of 14 SNe with spectro- 
scopic identifications or spectroscopic redshift determi- 
nations of their host galaxies, 8 objects are classified as 
case 1 (significant offsets and SN-like light curves) and 
5 objects are classified as case 3 (small offsets and SN- 
like light curves). Only 1 object are classified as case 4 
(significant offsets and non-SN-like light curves) proba- 
bly because the variable component of this object is the 
faintest among these 14 SNe (iyari — 25.3mag). Through 
our criterion for offsets of variability, 36% (= 5/14) of 
SNe are recognized as central variability, although spec- 
troscopic observations for SNe should be biased to off- 
set variability to avoid AGN. If we use this ratio, about 
88 (= 158 X 5/9) of 160 objects in case 3 are considered to 
be SNe. On the other hand, among 644 variable objects 
in these fields whose light curves covering longer than 
three years, the number of X-ray sources is 128. Out of 
128 optically variable X-ray sources - which we believe 
are AGN - 89 objects are in case 2, 37 objects are in case 
3, 1 object is in case 1, and 1 object is in case 4. Thus, 
30% (= (37 -f- 1)/128) of AGN have SN-hke fight curves. 
If we assume that this same fraction of all AGN have SN- 
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like light curves, then about 96 (= 228 x (37+l)/(89 + l)) 
of 160 objects in case 3 are considered to be AGN. The 
expected total number of these SNe and AGN in case 3 is 
184 (= 88 + 96) and reasonably consistent with the real 
number (160) in case 3. Hence, we made number densi- 
ties of SNe and AGN by scaling those for well-classified 
samples (case 1 for SNe and case 2 for AGN) by factors 
of 14/8 for SNe and 128/89 for AGN, as shown in Fig- 
ure [Tj] for the host object magnitudes and Figure [T51 for 
the variable component magnitudes. We discuss these 
figures in ^ We obtained these number densities of 489 
and 579 objects deg~^ for SNe and AGN, respectively. 

Nature of objects in case 4 is mysterious. We have 
only two spectroscopic data for objects in case 4 and 
both of the two spectra indicate that they are normal 
galaxies at z ~ 0.5. Their offsets from the centers of 
the host objects are significant. The light curves have 
marginal parameters around the threshold between those 
of SN-like light curves and non-SN-like light curves, and 
they may be misclassified as non-SN-like light curves due 
to measurement errors. On the other hand, only one 
object in case 4 is detected in X-ray. This object and 
some others have offsets from the centers of their host 
galaxies just above the threshold (1.2 pixel). Another 
possible reason why some objects are classified as case 
4 is misidentifications of the host objects due to their 
faintness. Such objects might be also transient objects. 
The number of objects in case 4 is not large and we do 
not include them for discussions of number densities. 

5. DETECTION COMPLETENESS OF OPTICALLY 
VARIABLE OBJECTS 

We intensively studied two types of completeness in 
this paper. One is completeness of variability detection 
itself, which should be functions of PSF sizes and back- 
ground noise of the images, and should be independent 
of properties of object variability. Another completeness 
depends on behaviors of object variability and observa- 
tional time samplings. It is not easy to estimate this 
type of completeness because properties of variability are 
complicated and differ from object to object. 

5.1. Variability Detection Completeness 

The first type of completeness is simply estimated by 
locating artificial variable objects (point sources with 
PSFs of the same size as those in the real images) ran- 
domly in the stacked images using the artdata task in 
IRAF, and detecting them in the same manner as for 
the real images. Three examples are shown in Figure 
m We can almost completely detect variable objects 
whose flux differences in magnitude unit are brighter 
than ~ 25.0 — 25.4 mag. The cut-offs at the bright 
ends around i'^^^^ = 20.0 — 20.5 mag are caused by mask- 
ing bright objects to avoid false detections of variabil- 
ity. The detection completeness decreases down to zero 
at i^j^j,; ~ 26.0 — 26.3 mag. The shapes of completeness 
curves are similar for all the subtracted images, although 
there are slight offsets in magnitude axis due to differ- 
ences of limiting magnitudes of the stacked images. The 
completeness is almost uniform in the whole regions of 
the subtracted images where we investigate object vari- 
ability. 
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Fig. 8. — Detection completeness of variable objects as a func- 
tion of variable component magnitude i'^^^^ calculated by adding 
artificial objects in the images and detecting them in the same 
manner as for the real images. Three cases are shown in the fig- 
ure. The limiting magnitudes of the images before subtractions are 
given between the top parentheses and drawn in dot-dashed lines. 
FWHM of PSF are also described between the bottom parenthe- 
ses. Cut-offs at bright magnitudes are caused by masking bright 
objects in order that our detection are not affected by saturation 
or non-linearity of the CCD pixel. 

5.2. Detection Completeness for Variable Stars, SNe, 
and ACN 

Variability detection completeness depends on shapes 
of light curves of variable objects and observational time 
samplings. We want to know detection completeness for 
each class of variable objects. It is difficult to do so 
for objects showing burst-like variability such as dwarf 
stars, SNe, and AGN, while we can easily simulate light 
curves of pulsating variable stars showing periodic vari- 
ability. AGN optical variability has been often charac- 
terized by the structure function (iKawaguchi et al.ill998l: 
Hawkinsl I2OOI Ide Vries eFall l2003t iVanden Berk et all 
200 llde Vries et al.l[2005l : ISesar et al.l[200l . and we can 
simulate AGN light curves to estimate the detection com- 
pleteness. Hence, we consider two types of simulated 
light curves; periodic variability for pulsating variable 
stars and variability characterized by the structure func- 
tions for AGN and estimate completeness. We also evalu- 
ate variability detection efficiency using real light curves 
and the X-ray source sample for AGN. 

In our time samplings and depths, SNe la up to z ~ 1.4 
can be detected in efficiency of a few tens of percent. 
Detection efficiency of SNe will be examined in detail in 
our following paper on SN rate. 

5.2.1. Periodic Variability 

Detection completeness for pulsating variable stars 
whose variability are periodic is easily calculated. Simu- 
lated light curves are parameterized by magnitude ampli- 
tudes A = 0.0 - 1.0 mag, periods T = 0.01 - 20000 days, 
and averaged magnitudes mg = 21 — 26 mag We have 
roughly two types of time baselines of observations, and 
then we show the completeness for the two fields, SXDF- 
C with observations from 2002 to 2005, and SXDF-N 
with those from 2002 to 2003, in Figure O There are 
slight differences between the detection completeness for 

■^^ Light curves are expressed as m{t) = ^ sin(27rt/T) + mo. 
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the two fields, but, both the completeness are almost the 
same on the whole. In all the five fields, some dark spikes 
indicating low sensitivity to those periods can be seen be- 
cause our observations have been carried out only in fall. 
The completeness is a strong function of magnitudes and 
amplitudes. We can detect only variable stars with large 
amplitudes {A > 0.4 mag) for faint stars of i' > 24 mag. 
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Fig. 9. — Detection completeness for variable objects showing 
periodic variability with periods T, amplitude A, and averaged 
magnitude mo in SXDS-C (10 epochs from 2002 to 2005, upper 6 
panels) and SXDS-N (8 epochs from 2002 to 2003, lower 6 panels). 



5.2.2. AGN Variability 

Unlike pulsating stars, AGN usually show burst-like 
variability aperiodically. Behaviors of AGN variability 
depends on rest-frame wavelength. Generally, AGN vari- 
ability are larger in shorter wavelength. Considering the 
time dilation of cosmological redshift, detection efficiency 
of AGN by optical variability in a certain broad band, in 
z'-band in this study, is not clear. To the zeroth ap- 
proximation, effects of wavelength dependence and time 
dilation are cancelled each other. 

Quasar variability has been often described in the form 
of the structure function. 



SF{At)= {m{t + At) - m{t)f /N{At) 



(1) 



where N{At) is the number of objects with time inter- 
vals of At. As far as we focus on time scales of months to 
several years (not several decades), a power-law form of 
the structure functions, SF{At) = {At/ro)^ , parameter- 
ized by a characteristic time scale of variability, tq, and 
a power-law slope, 7, is well fitted to observational data. 
The structure functions begin to show a turnover around 
the rest-frame time lag of a few years and this power- 
law approximation slightly overestirnates the variabil- 
itv (llvezic et al.ll2004U Ide Vries et "all 120051 : [Sesar et all 



[200 6[) . However, our observational data span over only 
three years at longest in observed frame and this simpli- 
fication is good enough for our estimation of AGN detec- 
tion completeness. 
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Fig. 10. — Detection completeness for variable objects whose 
variability are characterized by the structure functions in SXDS-C 
(10 epochs from 2002 to 2005, upper 6 panels) and SXDS-N (8 
epochs from 2002 to 2003, lower 6 panels). Typical values of 7 and 
TO for quasars are 0.2 and 6 X 10"" days, respectively, these values 
are plotted as stars. 



The slope of structure function provi des us impor- 
tant keys to the origi n of their variabilit y fHughes et al 
1992; HawkinSi 119961 : iKawaguchi et al. | [19 98: Hawkii^ 
20021) : accretion disk instability (iReesI 119841) . bursts of 
SN explosions (|Terlevich et aLl Il992f ). and microlens- 
ing ()Hawkinsll2002| ). Therefore, the slope of structure 



function has been investigated in many previous stud- 
ies. Moni toring observations of Paloma r- Green (PG) 
quasars ( Giveon et al.lll999l : lHawkinsl l2002) and variabil- 
ity studies of enormous number of SDSS quasars be- 
tween SDSS imaging data and older plate imaging data 
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(|de Vries et alJ [20031 [2001 iSesar et al.l[200l or SPSS 
spectrophotometric data (jVanden Berk et alJ [2004D set 
constraints on characteristic behaviors of variability. 
They obtained slopes of 7 ^ 0.2. We simulated light 
curves satisfying the structure functions with 5 x lO'' < 
To < 1 X 10^ days and 0.05 < 7 < 0.6 for averaged 
apparent magnitudes uiq — 21 — 26 mag for AGN com- 
ponents, and calculated detection completeness in our 
observational time samplings. We show results in Figure 
[TU|for two cases (4-year baseline for SXDF-C, -S, and -E 
and 2-year baseline for SXDF-N and -W) done for peri- 
odic variability in §5.2.11 Values of the slope obtained in 
previous studies indicate that main origin of AGN opti- 
cal variability is accretion disk instability or microlensing 
and the typical values of tq and 7 are ^ 6 x 10^ days and 
~ 0.2 in observed frame, respectively, which are indi- 
cated as stars in Figure These values were derived 
from observations in bands bluer than i'-band which is 
used in this work. In i'-band, it is expected that 7 stays 
roughly constant and tq become slightly larger. Deter- 
mination errors of tq and 7 are not small here and we 
infer that the completeness for our data is not less than 
~ 50% for mo < 24 mag. 

The detection completeness for quasars were also 
estimated using observational light curves of 42 PG 
quasars in B and j^ - bands over seven years obtained by 
iGiveon et al.l (|1999( l. The photometric points were well 
sampled and their observational time baselines are much 
longer than ours. They lacked time samplings in time 
scales of days and we interpolated the light curves. Their 
quasars are at relatively low redshift (0.061 < z < 0.371) 
with luminosity of —26.3 < Mb < —21.7 mag. The well- 
sampled light curves are very useful for calculating the 
completeness in our survey. In our estimations of the de- 
tection completeness, we assume the photometric errors 
of our surveys considering the time dilation and depen- 
dence of var iability on rest-frarn e wav elength which was 
indicated in IVanden Berk et al.l (|2004i) using the SDSS 
quasars. Variability depends on rest-frame wavelength 
that quasars are about twice as variable at 1000 A as at 
6OOOA as shown in their Equation 11 ^'^ and Figure 13. 
We obtained detection completeness curves as a function 
of observed i'-band magnitude as shown in Figure 1111 
From top left to bottom right, the redshifts are z = 0, 
1, 2, 3, 4, and 5. The results for the two fields, SXDF- 
C from 2002 to 2005 and SXDF-N from 2002 to 2003, 
are plotted. There are some differences between these 
cases and the longer baselines give us higher complete- 
ness. In this simulation, our detection completeness for 
AGN is down to zero at ~ 24.5 — 25.0 mag. Redshift 
dependence of the completeness are very small due to 
cancellation of time dilation and dependence of variabil- 
ity on rest-frame wavelength. We note that dependence 
of variability on AGN luminosity are not considered. In 
our detection limit, we can observe Seyfert-class AGN at 
low redshift and our estimates of completeness here may 
be lower limits for them because of the anti-correlation 
between AGN luminosity and variability. 

Optical variabihty of 172 X-ray sources are detected, 
while we have 481 X-ray sources in our variability survey 
field. The fraction of objects showing optical variability 
among X-ray sources is 36 ± 2%. All of the X-ray de- 

23 variability i;{A)[mag] = 0.616 exp(-A/988A) + 0.164. 
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Fig. 11. — Detection completeness for 42 P G quasars us- 
ing re al light curves over seven years obtained by IGiveon et al.1 
II1999I ). In the cases of objects at 2 = 0, 1, 2, 3, 4, and 5, 
we calculated the detection completeness considering cosmological 
time dilation and v ariability dependence on rest-frame wavelength 
llVand en Berk et al.ii2004 ') . We show the results for time samplings 
in SXDF-C in open circles and for those in SXDF-N in filled trian- 
gles. The fitted lines are l/[l-|-axexp{(mag — 6)/c}, parameterized 
by a, b, and c. 



tected variable objects are brighter than i' = 24.4 mag. 
Then, we also limit X-ray sources to those with i' < 24.4 
mag, the number is 334. Simply assuming that all the 
X-ray sources can show large variability enough to be 
detected, the detection completeness for X-ray sources 
is 51 ± 3%. Vanden B erk ct al. (2004) found that X- 
ray detected quasars show larger variability than quasars 
which are not detected in X-ray. This tendency is true 
for our sample; X-ray brighter sources show larger op- 
tical variability than X-ray fainter sources in a certain 
optical magnitude. Variable objects not detected in X- 
ray should have lower completeness than variable objects 
detected in X-ray. Some of the X-ray sources are type-2 
populations whose optical variability are more difficult 
to be detected given the unified scheme of AGN and de- 
tection completenss of type-1 AGN can be higher than 
these estimates. 

These three estimations of detection completeness for 
AGN are roughly consistent with each other and strongly 
depend on apparent magnitudes as well as that for pul- 
sating variability. Inferred completeness is ~ 1 at i' ~ 21 
mag and slowly decreases down to zero at = 24 — 25 
mag with uncertainty of a factor of a few. 

6. OVERALL SAMPLE 

We found 1040 variable objects including possible tran- 
sient objects in i'-band in the SXDF down to variable 
component magnitudes of i'^^^.^ ^25.5 mag. In the top 
left panel of Figure [121 the number densities N of all 
the detected variable objects are shown as a function of 
time interval At of observations when we compare only 
two separated images. We also draw a line fitted in the 
form of N{At) = a[l — exp{—b{Aty}] Number densities 
as a function of At for SNe are expected to bend over at 
a certain time interval because variability time scales of 
SNe are roughly common to each other, a few months, in 
rest-frame. Then, we introduce such function forms for 
all the kinds of variable objects. Detection limits for vari- 
able objects are slightly different for each combination of 



Optically Faint Variable Object Survey in the SXDF 



13 



observational epochs and it can affect obtained number 
densities. Then, we also show number densities of vari- 
able objects with variable components brighter than 25.0 
mag in the right panels. Monotonic increases of number 
densities in observed frame indicate that objects show- 
ing variability in time scales of years are dominant in 
this magnitude range (and in the SXDF direction; the 
Galactic halo). The number densities for variable stars 
are shown in the second row, SNe in the third row, and 
AGN in the bottom row after classifying variable objects 
in IjH Variable objects for calculating number densities 
used in this figure are 644 objects with light curves in 
long baselines from 2002 to 2005 enough to discriminate 
SNe from AGN, and are within the IRAC field. These 
SN and AGN number densities are corrected for case 
3 objects. Typical time scales of variability are from 
days to months for variable stars and SNe while AGN 
show variability in time scales of months to years (e.g., 
iVanden Berk et ani2004[) . The flat distributions for vari- 
able stars indicate that stars showing variability in short 
time scales are mainly included in the sample. The SN 
number densities arrive at plateau in time scales of about 
a few months as expected from time scales of SN light 
curves. On the other hand, significant increases towards 
time scales of years are seen for AGN. Fitted lines for 
AGN in the same form of function as for SNe show pos- 
sible tur novers at ~ 1000 days. Results of SDSS 
quasars (jlvezic et al.ll2003l ) indicated turnovers of struc- 
ture functions at At ~ 1000 days in rest-frame and these 
possible turnovers in our results may be real. These re- 
sults indicate that long time baselines of years make the 
completeness significantly higher only for AGN. 
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Fig. 12. — Number densities of variable objects as a function 
of time intervals of At in i'-band imaging observations. We plot 
them for objects showing variability above bcri^ in the left column, 
objects with variable component magnitude i'^^^.^ brighter than 25.0 
mag in the right column. These points are fitted in the form of 
N = a X exp{—b{At)'^} (solid lines). The fitted line in the top left 
panel is shown in gray lines in every panel as a reference. 



We show the number densities of variable objects as 
a function of variable component magnitude i^^^.; in the 
left column of Figure [131 Number densities after classi- 



fying variable objects into variable stars, SNe, and AGN 
in three time scales. At < 10 days in the second column, 
10 < At < 50 days in the third column, 50 < At < 200 
days in the fourth column, and At > 200 days in the 
right column, are also plotted. All the number densities 
drop around i'^^^.^ ~ 25 mag and our variability detections 
are reasonably consistent with simulated completeness as 
shown in Figure [S] In the bottom row of Figure [131 the 
number densities of transient objects are plotted. By def- 
inition of transient objects (see [21) i variable component 
magnitudes of them correspond to total magnitudes in 
their brighter phases. The transient object sample can 
include not only objects showing transient phenomena 
such as flare-ups of faint dwarf stars but also less slowly 
moving Kuiper belt objects than uj ^ 1" hour~^ (a typi- 
cal value of seeing size per exposure time for each stacked 
image), corresponding to a semimajor axis of > 100 AU. 
These estimates of numbers of variable objects provide us 
the expected numbers of variable objects contaminated 
into interested samples u sing non -simultaneo us observa- 
tional data. For example, live et al. (2006 ) and lOta et al.l 
(|2007( ) investigated the possibility of variable object con- 
taminations into thei r LAE sample at z ^ 7. in the Sub- 
aru Deep Field fSDF^ lKashikawa et al.ll200"^ . which were 
obtained comparing broad-band imaging data taken be- 
fore 2004 with narrow-band NB973 imaging data in 2005. 
From the bottom row of Figure [T3l a few or less than one 
objects in the SDF (~ 0.3 deg^) can be just transient ob- 
jects and misclassified as narrow-band excess objects in 
their sample. Their plausible candidates with narrow- 
band excess are two. One of them were spectroscopically 
identified and turned out to be a real LAE at z = 6.94. 
Spectroscopic identification for another candidate should 
be done. Whether this candidate is a real LAE or not, 
the number of narrow-band excess objects they found is 
consistent with expected number from the statistics of 
number densities of transient objects. 

The fractions of variable objects to the overall objects 
in the SXDF are shown in Figure [141 Magnitudes used 
in this figure are the total magnitudes of host objects. 
About 5% of objects at i' ^ 21 mag show optical vari- 
ability, and the fractions rapidly decrease towards fainter 
magnitudes. These declines are caused by detection limit 
for variability. Large variability relative to the host ob- 
jects are necessary to be detected for faint objects. Al- 
most all of objects in the SXDF are galaxies, not stars 
(the fraction of stars is ^ 4%). Then, the fraction of 
variable AGN is ~ 3% around i' ^ 21 mag, wh e re the 
detection completeness is ~ 1. ISaraiedini et al.l (|2006f ) 
found that 2.6% of galaxies have variable nuclei down to 
^nuc < 27.0 mag without completeness corrections using 
two-epoch observations separat ed by seven y ears. Their 
fraction is consistent with ours. ICohen et all (2006) also 
used four-epoch ACS imaging data in time baselines of 
three months down to Vtotai < 28 mag and detected vari- 
ability of ~ 1% o f galaxies. This small percentage in 
ICohen et al.l ()2006^ is consistent with ours if we consider 
their short time baselines because most of AGN vary in 
brightness in longer time scales of years as shown in Fig- 
ure [HI 

7. VARIABLE STARS 

The sample of variable stars used in this section in- 
cludes 153 objects which were classified as 1) reliable 
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Fig. 13. — Number densities of variable objects as a function 
of variable component magnitude i^j^^j. We plot those for all the 
variable objects in the top row, variable stars in the second row, 
SNe in the third row, and AGN in the bottom row, respectively. 
Number densities in each time scale, Ai < 10 days in the second 
column, 10 < At < 50 days in the third column, 50 < At < 200 
days in the fourth column, and At > 200 days in the right column, 
as well as in the case using all the observational data in the left 
column, are also plotted. Number densities of possible transient 
objects are shown in the bottom row. 
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Fig. 14. — Fractions of all the variable objects to the whole 
objects in the SXDF as functions of i'-band magnitude, B — V, 
and R — i' colors in the top row. Fractions after classifying variable 
objects are also plotted in the second, third, and bottom rows. 



stars, 2) probable stars, and 3) possible stars in §4.11 
The top panels of Figure [15] show color-magnitude di- 
agrams of i' versus i? — i' in the left panel and V versus 
i? — in the right panel for variable stars (black cir- 
cles) and non- variable stars in the SXDF (gray dots). 
These non-variable stars are selected through the same 
criteria as for the variable stars. Number counts and 
fractions of variable stars are also plotted in the rows be- 
low. Figure [16] shows color distributions and fractions of 
variable stars for R — i' and B — V . The number count 
in i'-band has a relatively steep peak at i' ~ 21.5 mag 
and decreases towards fainter magnitudes. This cut-off 
is partly caused by the variability detection limit. The 
dot-dashed lines indicated in the left columns of Figure 
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Fig. 15. — Top panels show distributions of variable objects 
(black circles) in color-magnitude diagrams, V versus B — V in the 
left panel and i' versus R — i' in the right panel. Non- variable stars 
are also shown in gray dots. The criterion for excluding contami- 
nations of galaxies, B — V > 0.08 y.V — 1.59, is drawn in a dotted 
line in the V versus B — V diagram. Number counts and fractions 
of variable objects for i' and V magnitudes are shown in the second 
and third rows, respectively. 
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Fig. 16. — Histograms of numbers of variable objects are shown 
in the upper panels. The fractions are also plotted in the lower 
panels. The left column is for R — i' color and the right column is 
for B-V color. 



[T5I are total magnitudes of objects with variable com- 
ponents of the detection limit, i'^^^.-^ ~ 25.5 mag, in the 
cases of variability amplitudes of 0.05, 0.1, and 0.3 mag 
from left to right. Thus, limiting magnitudes can be as 
shallow as ^ 22.2 mag for variable objects with low vari- 
ability amplitudes of ^ 0.05 mag. Poor time samplings of 
our observations prevent us from determining amplitudes 
of variability in cases of periodic variability. Assuming 
that amplitudes roughly equal to magnitude differences 
of objects between their maxima and minima in our sam- 
plings, amplitudes of almost all of variable stars are less 
than 0.1 mag as shown in Figure [T71 These variable stars 
with low amplitudes can be detected above ~ 22.2 mag. 
Variable stars with larger amplitudes can be detected 
for fainter stars. The top right panel of Figure [T7I clear Iv 
indicates this tendency of variability selection effects. 

Bimodal distributions are clearly seen in the figures of 
color-magnitude diagrams (FigurefTSj) and color distribu- 
tions (Figure fTB)) : there are two populations of blue bright 
(y ~ 22 mag) variable stars and red faint {V ~ 23.5 
mag) variable stars. From comparison with the Besangon 
model prediction, stars in the upper sequence belong to 
the thick disk populations rather than thin disk because 
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stellar sequence of the thin disk is expected to be redder 
than the observed sequence hy B — V of ^ 0.3. The red 
variable stars are considered to be dwarf stars showing 
flares or pulsating giant stars. On the other hand, stars 
in the lower sequence belong to the halo populations. 
The blue variable stars are considered to be pulsating 
variable stars such as S Scuti, RR Lyrae, and 7 Doradus 
stars, within the instability strip in color-magnitude di- 
agrams. Eclipsing binaries may be also included. 





1 


1 2 20 22 24 26 28 
R—i' [mog] i'— band [mag] 







1 2 20 22 24 26 28 

B-V [mag] K-bond [mag] 



Fig. 17. — Variability amplitudes of variable stars as function of 
apparent magnitude (i'-band and V-band) and colors {R — i' and 
B — V). These amplitudes are defined as magnitude differences 
between maxima and minima in our time samplings and give lower 
limits of real amplitudes. 

Ilvezic et"aLl (|200C1[) selected RR Lyrae candidates us- 
ing SDSS multi-epoch wide-field imaging data in short 
time scales of days, and found a cut-off magnitude at 
r* ~ 21 mag. They indicated that this cut-off magni- 
tude corresponds to a distance of ^ 65 kpc from the 
Galactic center, which might be the edge of the Galac- 
tic halo. Our time samplings are too poor to determine 
their periods of variability as well as real amplitudes. 
Maxima of magnitude differences in our time samplings 
shown in Figure [T7] are expected to be smaller than real 
amplitudes by a factor of 1.2 from simulated light 
curves. Our typical exposure time of one hour, which 
is not much lo nger than va. riability periods of RR Lyrae 
(0.3 — 0.6 davs. lVivas et al.|[2004 ). also makes these mag- 
nitude differences small. Moreover, variability in ?'-band 
are smaller than those in l^-band for RR Lyrae (0.5 — 1.0 
mag in V^-band, Vivas et al. 20(33) b y a factor by ~ 1. 5 
(interpolation of values in Table 7 in lLiu fc Janesl[T990D . 
In total, the maxima of magnitude differences in i'-band 
are expected to be less than 0.2 — 0.5 mag for RR Lyrae 
in our time samplings. Given the variability amplitudes 
of blue variable stars with B — V 0.3 in our sample, 
there are one or a few RR Lyrae candidates with am- 
plitudes as large as those of RR Lyrae (the bottom left 
panel of Figure [17]). Variability amplitude of the spectro- 
scopically identified blue variable star with B — V = 0.44 
shown in the bottom panel of Figure [3] is small and can 
not be a candidate of RR Lyrae. If these RR Lyrae can- 
didates in our sample are really RR Lyrae, the distances 
from the Galactic center, which are calculated by appar- 
ent y-band magnitude, are larger than ^ 150 kpc. The 
inferred number densities of RR Lyrae are ^ 10~^ kpc~'^ 
at a distance of 150 — 250 kpc from the Galactic center 
towards t he halo, and consi stent with extrapo l ations of 
results bv llvezic et all (|2000D and lVivas fc ZinnI (|2006D in 



spite of our poor statistics. They are candidates of the 
most distant Galactic stars known so far. Other types of 
blue variable stars in the halo population than RR Lyrae 
are considered to be included in the sample such as Popu- 
lation II 6 Scuti stars, 7 Doradus stars, eclipsing binaries, 
and so on. In order to examine nature of these faint blue 
variable stars, dense monitoring observations for deter- 
minations of pulsating periods and amplitudes, and/or 
follow-up spectroscopic observations are necessary. 

Two of the red variable stars with B — V = 1.36 and 
1.34 were spectroscopically observed and identified as M 
dwarf stars (Figure [3|). The B — V colors of the red 
faint variable stars are ~ 0.8 — 2.0, which are those of 
K or M type stars. The F-band absolute magnitudes of 
these stars are widely spread, for example. My = 8 — 
12 mag for M stars with colors of 1.3 < B — < 1.7 
in the case of ma in sequence from the Hipparcos data 
(|Koen et al.ll200^ . Assuming that almost all of these 
red variable stars are dwarf stars, not red giant stars, 
the inferred distances to the stars are ^ 1.5 — 4 kpc from 
the Galact ic plane, where the thick disk population is 
dominant fiChen et al.ll200l! ) . The number density of the 
variable dwarf stars is ~ 2 x 10'^ kpc~^. Magnitude and 
color distributions of stars from the model predictions are 
consistent with those stars in the SXDF selected through 
our criteria and it might indicate that a few percent of 
the whole dwarf stars show optical bursts in our time 
samplings as shown in the bottom rows of Figure 1161 

Fractions of variable stars are one of interesting results. 
Fractions of variable stars as functions of colors oi R~i' 
and B — V are nearly fl a t at a few percents as shown in 
Figure [161 iTonrv et al.l ()2005 ) investigated fractions of 
variable stars as a function of quartile variability down 
to variability of 0.015 mag. They examined variability of 
stars in sequential 14-day observations down to i? ~ 20 
mag at a superb photometric precision of 0.002 mag. 
Their result indicates that a fraction of 0.016 mag/x of 
stars show variability above quartile variability of x mag. 
Our detection threshold for variability of point sources 
corresponds to about 0.05 — 0.10 mag and expected frac- 
tions are ~ 0.02 — 0.03 using this equation. The fraction 
obtained by our survey for blue stars is ~ 0.03 while that 
for red stars is ~ 0.05 on average, and are slightly larger 
than the expected fraction. This discrepancy might be 
derived from some reasons des cribed below . A y oung 
open cluster NGC2301, which iTonrv et all (j2005f ) tar- 
geted, is only 146 Myr old, where there should not be 
many variable populations. On the other hand, our 
pointings are towards the halo, where main components 
observed are old populations. Another reason can be the 
difference of time samplings; their observations are 14 
consecutive days while our observations are sparse over 
years. Time scales which can be examined by each study 
are clearly different, which can cause these differences. 
FSVS studies recently investigated fractions of faint vari- 
able stars down to V ~ 24 mag (jMorales-Rueda et al.l 
l2006t iHuber 'etelW m)^ . Their survey depths and direc- 
tions are similar to ours , and their results for fractions of 
variable stars, 5 — 8% (jHuber et al.ll200"6l ) are also con- 
sistent with ours. 

8. SUMMARY 

We investigated optical variability of faint objects 
down to i^jjj.; 25.5 mag for variable components over 
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0.918 deg2 in the SXDF. Multi-epoch (8-10 epochs over 
2 — 4 years) imaging data obtained with Suprime-Cam 
on Subaru 8.2-m telescope provided us the first statisti- 
cal sample consisting of 1040 optically variable objects by 
image subtraction for all the combinations of images at 
different epochs. We classified those variable objects into 
variable stars, SNe, and AGN, based on the optical mor- 
phologies, magnitudes, colors, optical-mid-infrared col- 
ors of the host objects, spatial offsets of the variable 
components from the host objects, and light curves. Al- 
though not all the variable objects were classified because 
of short time baselines of observations in the two fields, 
our classification is consistent with spectroscopic results 
and X-ray detections. 

We examined detection completeness for periodic vari- 
ability and AGN variability. The completeness strongly 
depends on apparent magnitude. The completeness for 
AGN is '--^ 1 at i' ~ 21 mag and deceases down to zero 
at i' ^ 24.5 mag. Rcdshift dependence of the complete- 
ness calculated using light curves of PG quasar is small 
due to cancellation of time dilation and anti-correlation 
between rest-frame wavelength and variability. Among 
X-ray sources in the field, 36 ± 2% (51 ± 3% for the 
bright sample with i' < 24.4 mag) show optical vari- 
ability. Variability detections of X-ray sources also show 
similar dependence on apparent magnitude to that from 
light curves. 

Number densities of variable objects for the whole sam- 
ple, variable stars, SNe, and AGN as functions time in- 
terval At and variable component magnitude i'^^^.^ were 
obtained. About 5% of all the objects show variability at 
i' = 21 — 22 mag including host components although de- 
creasing fractions towards fainter magnitude are caused 
by the detection limit for variable components. Number 
density of variable stars as a function of time interval At 
is almost flat indicating that time scales of variability of 
these stars are short. Number density of SNe arrives at 
platau in time scales of a few months and that of AGN 
increases even in time scales of years. These results are 
consistent with expectations from typical time scales of 
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their variability. Total number densities of variable stars, 
SNe, and AGN are 120, 489, and 579 objects deg-^, re- 
spectively. 

Variable stars show bimodal distributions in the color- 
magnitude diagrams. This indicates that these variable 
stars consist of blue bright {V ^ 22 mag) variable stars 
of the halo population and red faint {V ~ 23.5 mag) vari- 
able stars of the disk population. We selected a few can- 
didates of RR Lyrae considering their large magnitude 
differences between maxima and minima and blue B — V 
colors. The number density is ^ 10^^ kpc^'^ at a distance 
of > 150 kpc from the Galactic center, which is consis- 
tent with extrapolations of previous results. Follow-up 
observations to determine the amplitudes and periods of 
variability might show that these candidates are at the 
outermost region of the Galactic halo. 

Our statistical sample of optically variable objects pro- 
vides us unique opportunity for the studies such as AGN 
properties and SN rate, which will be topics in our fol- 
lowing papers. There are planned large surveys such as 
Panoramic Surv ey Teles cope and Rapid Response Sys- 
tem fPan-Starrs HKaiseT ct al. 200^, Large Synoptic Sur- 
vey Telescope fLSST: [Tvsqn 20021). and SuperN ova Ac- 
celeration Probe fSNAP: lAldering et al.] I2002D . This 
work also provides basic information for such future wide 
and deep variability surveys. 
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TABLE 1 

Imaging Data Other than Optical Variability 



Wavelength 


Telescope /Instrument 


Band 


Detection Limit 


Area'' [deg^] 


Variable'' 


Detection'^ 


optical 


Subaru/Suprime-Cam'^ 


B,V,Rc,i',z 


28.2,27.2,27.6,27.5,26.5 mag<= 


0.918 


1040 


1040 


X-ray 


XMM-Newton/EPIC 


2.0-10.0 keV 


3 X 10~^^ erg s"^ cm~^ 


0.808 


936 


165 


X-ray 


XMM-Newton/EPIC 


0.5-2.0 keV 


1 X lO"'^^ erg s~^ cm~^ 


0.808 


936 


91 


mid-infrared 


Spitzer/IRAC 


3.6^m-band 


22.0 mag'^ 


0.889 


1028 


995 



^ Area overlapping with the variability survey field. ^ Number of optically variable objects within the arca.*^ Number of optically 
variable objects detected in the bands within the area. Derived from preliminary version of the SXDF catalogs in IFurusawa et al.l 
(120071) .^ Detection limits for optical and mid-infrared imaging arc measured in AB magnitude. 
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TABLE 2 

Summary of Subaru Suprime-Cam i'-band Imaging Observations 

FOR A'ARIAniLIT^' DfTFCTIONS 



field epoch date (UT)* At^ texp [sec] mum'^ 



SXDF-C 


1 


02/09/29,30 


0.0 


2700 


0.54 


26.19 


SXDF-C 


2 


02/11/01 


32.6 


1860 


0.92 


25.76 


SXDF-C 


3 


02/11/02 


33.5 


1800 


0.68 


25.85 


SXDF-C 


4 


02/11/05 


36.7 


2400 


0.70 


26.11 


SXDF-C 


5 


02/11/09 


40.5 


2460 


0.60 


25.77 


SXDF-C 


6 


02/11/27,29 


59.8 


4200 


0.72 


26.38 


SXDF-C 


7 


02/12/07 


68.4 


3000 


0.78 


26.32 


SXDF-C 


8 


03/10/20 


385.7 


5760 


1.14 


26.53 


SXDF-C 


9 


03/10/21 


386.5 


7500 


0.58 


26.71 


SXDF-C 


10 


05/09/28 


1094.6 


3600 


1.00 


26.04 


SXDF-N 


1 


02/09/29,30 


0.0 


3300 


0.56 


26.26 


SXDF-N 


2 


02/11/01 


32.4 


2640 


0.96 


25.88 


SXDF-N 


3 


02/11/02 


33.3 


1800 


0.68 


25.86 


SXDF-N 


4 


02/11/09 


40.3 


2100 


0.64 


25.78 


SXDF-N 


5 


02/11/29 


60.3 


3300 


0.74 


26.27 


SXDF-N 


6 


03/09/22 


357.5 


4264 


0.60 


26.37 


SXDF-N 


7 


03/10/02 


367.6 


1500 


0.70 


25.88 


SXDF-N 


8 


03/10/21 


386.5 


3000 


0.72 


26.14 


SXDF-S 


1 


02/09/29,30 


0.0 


3000 


0.52 


26.28 


SXDF-S 


2 


02/11/01 


32.5 


3600 


1.04 


25.91 


SXDF-S 


3 


02/11/02 


33.4 


1800 


0.70 


25.83 


SXDF-S 


4 


02/11/09 


40.4 


2580 


0.66 


25.60 


SXDF-S 


5 


02/11/29 


60.6 


1500 


0.82 


26.00 


SXDF-S 


6 


03/09/22 


357.6 


4500 


0.54 


26.45 


SXDF-S 


7 


03/10/02 


367.7 


2040 


0.68 


26.00 


SXDF-S 


8 


05/09/28 


1094.6 


3900 


0.96 


26.04 


SXDF-E 


1 


02/09/29,30 


0.0 


3000 


0.60 


26.25 


SXDF-E 


2 


02/11/01 


32.5 


3000 


1.04 


25.97 


SXDF-E 


3 


02/11/02 


33.4 


1800 


0.70 


25.83 


SXDF-E 


4 


02/11/09 


40.4 


2820 


0.66 


26.78 


SXDF-E 


5 


02/11/29 


60.6 


1800 


0.80 


26.03 


SXDF-E 


6 


02/12/07 


68.6 


1209 


1.54 


25.19 


SXDF-E 


7 


03/09/22 


357.6 


6000 


0.62 


26.58 


SXDF-E 


8 


03/10/02 


367.7 


1271 


0.68 


25.71 


SXDF-E 


9 


03/10/21 


386.7 


1400 


0.88 


25.51 


SXDF-E 


10 


05/09/28 


1094.6 


3600 


0.96 


26.11 


SXDF-W 


1 


02/09/29,30 


0.0 


2400 


0.54 


26.14 


SXDF-W 


2 


02/11/01 


32.7 


3000 


0.96 


25.96 


SXDF-W 


3 


02/11/02 


33.5 


1800 


0.66 


25.84 


SXDF-W 


4 


02/11/05 


36.8 


3060 


0.76 


25.98 


SXDF-W 


5 


02/11/09 


40.5 


2100 


0.64 


25.97 


SXDF-W 


6 


02/11/27,29 


59.8 


4200 


0.74 


26.39 


SXDF-W 


7 


02/12/07 


68.5 


6483 


1.04 


26.34 


SXDF-W 


8 


03/10/20,21 


386.5 


5460 


0.66 


26.46 



Observed date in yy/mm/dd. When the images were 
stacked together with those at different dates, both dates 
are included.*' Days from the first observation in each 
field.'' FWHM of PSF in stacked images.'' Limiting mag- 
nitude of 5a in 2'.'0 aperture. 



